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Abstract. - The magnetoresistance (MR) of CeColns is notably different from that expected 
for orbital MR due to the Lorentz force and described by Kohler's rule which holds in many 
conventional metals. We show that a pronounced crossover from negative to positive MR of 
CeColns that occurs at elevated temperatures is determined by the dependence of the effective 
mass M* (B, T) on both magnetic field B and temperature T. Thus, the crossover is regulated by 
the universal behavior of M* (B, T) observed in heavy-fermion metals. This behavior is exhibited 
by M* (B,T) when a strongly correlated electron system transits from the Landau Fermi liquid 
behavior induced by the application of magnetic field to the non-Fermi liquid behavior taking 
place at rising temperatures. Our calculations of MR are in good agreement with facts and reveal 
new scaling behavior of MR. 
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An explanation of the rich and striking behavior of 
strongly correlated electron system in heavy fermion (HF) 
metals is, as years before, among the main problems of 
the condensed matter physics. One of the most interest- 
ing and puzzling issues in the research of HF metals is the 
anomalous normal-state transport properties. HF metals 
show a number of distinctive transport properties, among 
which is the magnetoresistance (MR) of HF metals. Mea- 
surements of MR on CeColns [1,2] have shown that this is 
notably different from that expected for weak-field orbital 
MR and described by Kohler's rule which holds in many 
conventional metals, see e.g. [3]. MR of CeColns exhibits 
a crossover from negative to positive MR that occurs in 
fixed magnetic fields B with increasing temperature T, so 
that at the high fields and relatively low temperatures MR 
becomes negative [1,2]. 

This crossover is hard to explain within the conventional 
Fermi liquid theory for metals and in terms of Kondo sys- 
tems [4] and therefore it is assumed that the crossover can 
be attributed to some energy scales causing a change in 
character of spin fluctuations with increasing the applied 
magnetic field strength [1] . It is widely believed that such 
quantum fluctuations becoming sufficiently strong sup- 
press quasiparticles at a quantum phase transition and 
when the system in question transits from its Landau- 
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Fermi liquid (LFL) regime to non-Fermi liquid (NFL) be- 
havior [5-7]. 

On the other hand, even early measurements carried out 
on HF metals gave evidences in favor of the existence of 
quasiparticles. For example, the application of magnetic 
field B restores LFL regime of HF metals which in the 
absence of the field demonstrates NFL behavior. In that 
case the empirical Kadowaki- Woods ratio K is conserved, 
K = A(B)/~f%(B) oc A(B)/ X 2 {B) = const [8,9] where 
C/T = 7o, C is the heat capacity, \ is magnetic suscep- 
tibility and A(B) is the coefficient determining the tem- 
perature dependence of the resistivity p = p n + A(B)T 2 . 
Here po is the residual resistance. The observed conser- 
vation of K can be hardly accounted for within scenarios 
when quasiparticles are suppressed, for there is no rea- 
son to expect that jo{B), x(T), A(B) and other transport 
and thermodynamic quantities like the thermal expansion 
coefficient a(B) arc affected by the fluctuations or local- 
ization in a correlated fashion. 

Quasiparticles were observed in LFL regime in measure- 
ments of transport properties on CeColns [10]. While it 
is extensively accepted that the NFL behavior is deter- 
mined by the critical fluctuations, Kondo lattice [5-7] and 
multiple energy scales [11], therefore in that scenario the 
crossover region has to be formed by the fluctuations and 
scales rather then by quasiparticles. Analyzing thermody- 
namic quantities, it was shown that quasiparticles exist in 
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both the LFL and the crossover regimes when strongly cor- 
related Fermi systems such as HF metals [12-16] or two- 
dimensional 3 He [17] transit from its LFL to NFL behav- 
ior. Therefore it is of crucial importance to verify whether 
quasiparticlcs characterized by their effective mass M* 
still exist and determine the transport properties of HF 
metals in the crossover region. As we will see, measure- 
ments of MR in the crossover region can present indica- 
tive data on the availability of quasiparticles. Fortunately, 
such measurements of MR were carried out on CeCoIns 
when the system transit from the LFL to NFL behavior 
at elevated temperatures and fixed magnetic fields [1,2]. 

In this Letter we analyze MR of CeCoIns and show 
that the crossover from negative to positive MR that oc- 
curs at elevated temperatures and fixed magnetic fields 
can be well captured utilizing fermion condensation quan- 
tum phase transition (FCQPT) based on the quasiparti- 
cles paradigm [14,18-20]. We demonstrate that crossover 
is regulated by the universal behavior of the effective mass 
M*(B,T) observed in many heavy-fermion metals and is 
exhibited by M*(B,T) when HF metal transits from the 
LFL behavior induced by the application of magnetic field 
to NFL behavior taking place at rising temperatures. Our 
calculations of MR are in good agreement with facts and 
allow us to reveal new scaling behavior of MR. Thus, we 
show that the transport properties are mainly determined 
by quasiparticlcs rather then by the critical fluctuations, 
Kondo lattice and energy scales which are expected to ar- 
range the behavior in the transition region. 

To study universal low temperature features of HF met- 
als, we use the model of homogeneous heavy-fermion liq- 
uid with the effective mass M*(T, B,x), where the number 
density x = p F /3ir 2 and pp is the Fermi momentum [21]. 
This permits to avoid complications associated with the 
crystalline anisotropy of solids [15]. We first outline the 
case when at T — > the heavy-electron liquid behaves as 
LFL and is brought to the LFL side of FCQPT by tun- 
ing a control parameter like x. At elevated temperatures 
the system transits to the NFL state. The dependence 
M*(T,x) is governed by Landau equation [21] 
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where n(p, T, x) is the distribution function of quasiparti- 
cles and F(pf, p) Landau interaction amplitude. At T = 
0, eq. (1) reads [21] M*/M = 1/(1 - N F 1 (p f ,Pf)/3). 
Here Nq is the density of states of a free electron gas, 
F 1 (pf,pf) is the p-wave component of Landau interac- 
tion amplitude F. Taking into account that x = p F /3ir 2 , 
we rewrite the amplitude as F 1 (pf,Pf) — F 1 (x). When 
at some critical point x = xpe, F (x) achieves certain 
threshold value, the denominator tends to zero and the 
system undergoes FCQPT related to divergency of the ef- 
fective mass [14,18,20] 
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where M is the bare mass, eq. (2) is valid in both 3D and 
2D cases, while the values of factors A and B depend on 
the dimensionality. The approximate solution of eq. (1) 
is of the form [16] 
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where Ai > and A2 < are constants of order unity, 
P = TM*(T)/p% and /(0) - F 1 {x FC )- It follows from 
eq. (3) that the effective mass M* as a function of T 
and x reveals three different regimes at growing temper- 
ature. At the lowest temperatures we have LFL regime 
with M*(T,x) ~ M*(x) + aT 2 with a < since A : > 0. 
The effective mass as a function of T decays down to a 
minimum and after grows, reaching its maximum M^j at 
some temperature Tm(x) then subsequently diminishing 
as T -2 / 3 [12, 14]. Moreover, the closer is the number den- 
sity x to its threshold value xpe, the higher is the rate of 
the growth. The peak value grows also, but the max- 
imum temperature Tm lowers. Near this temperature the 
last "traces" of LFL regime disappear, manifesting them- 
selves in the divergence of above low-temperature series 
and substantial growth of M*(x). Numerical calculations 
based on eqs. (1) and (3) show that at rising temperatures 
T > T1/2 the linear term cx f3 gives the main contribution 
and leads to new regime when eq. (3) reads M/M* (T) oc (3 
yielding 

M*(T) oc T~ 1/2 . (4) 

We remark that eq. (4) ensures that at T > Ti/ 2 the 
resistivity behaves as p(T) oc T [14]. 
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Fig. 1: Schematic plot of the normalized effective mass = 
M*(T/T m )/Mm {Ml, is its maximal value at T = T M ) versus 
the normalized temperature Tjv = T/Tm- Several regions are 
shown. First goes the LFL regime (Mjv(TW) ~ const) at Tn <S 
1, then the transition regime (the hatched area) where M% 
reaches its maximum. At elevated temperatures T _2//3 regime 
occurs followed by T~ 1//2 behavior, see eq. (4). 

Near the critical point xfc, with M/M*(x — * xfc) — * 
0, the behavior of the effective mass changes dramatically 
because the first term on the right-hand side of eq. (3) 
vanishes, the second term becomes dominant, and the ef- 
fective mass is determined by the homogeneous version of 
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(3) as a function of T. As a result, the scale M/M* van- 
ishes and we get to scale M* in M M and T in Tm- These 
scales can be viewed as natural ones. The schematic plot 
of the normalized effective mass = M* /M M versus 
normalized temperature Tn = T '/Tm is reported fig. 1. In 
fig. 1 both T -2 / 3 and T -1 / 2 regimes are marked as NFL 
ones since the effective mass depends strongly on temper- 
ature. The temperature region T ~ Tm{x) signifies the 
crossover between the LFL regime with almost constant 
effective mass and NFL behavior, given by T~ 2 / 3 depen- 
dence. Thus temperatures T ~ Tm can be regarded as 
the crossover between LFL and NFL regimes. It turns 
out that M*(T,x) in the entire T < T" 1 / 2 range can 
be well approximated by a simple universal interpolating 
function [12,14,16]. The interpolation occurs between the 
LFL (M* oc T 2 ) and NFL (M* oc T" 1 / 2 , see cq. (4)) 
regimes thus describing the above crossover. Introducing 
the dimcnsionlcss variable y = TV = T '/Tm, we obtain the 
desired expression 
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Here M^(y) is the normalized effective mass, c\ and C2 
are parameters, obtained from the condition of best fit to 
experiment. To correct the behavior of M^(y) at rising 
temperatures we add a term to eq. (5) and obtain 
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where C3 is a parameter. The last term on the right hand 
side of eq. (6) makes satisfy eq. (4) at rising temper- 
atures T/T M > 2. 

At small magnetic fields B (that means that the Zeeman 
splitting is small), the effective mass does not depend on 
spin variable and B enters eq. (1) as B^b/T making 
Tm oc B\xb where fis is the Bohr magneton [14,16]. The 
application of magnetic field restores the LFL behavior, 
and at T = the effective mass depends on B as [12, 14] 
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where B c q is the critical magnetic field driving both HF 
metal to its magnetic field tuned QCP and corresponding 
Neel temperature toward T = 0. In some cases B c q = 0. 
For example, the HF metal CeRu2Si2 is characterized by 
B c q = and shows neither evidence of the magnetic or- 
dering or superconductivity nor the LFL behavior down 
to the lowest temperatures [22]. In our simple model B c q 
is taken as a parameter. We conclude that under the ap- 
plication of magnetic filed the variable 



y = T /Tm oc 



T 



H B (B-B c0 ) 
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remains the same and the normalized effective mass is 
again governed by eqs. (5) and (6) which arc the final 
result of our analytical calculations. We note that the 



obtained results are in agreement with numerical calcula- 
tions [12,14]. 

The normalized effective mass Mjy(y) can be ex- 
tracted from experiments on HF metals. For example, 
M*(T,B) oc C{T)/T oc S(T)/T oc xac{T), where S(T) 
is the entropy, C(T) is the specific heat and xac{T) is 
ac magnetic susceptibility. If the corresponding measure- 
ments are carried out at fixed magnetic field B (or at fixed 
both the concentration x and B) then, as it seen from fig. 
1, the effective mass reaches the maximum at some tem- 
perature Tm- Upon normalizing both the effective mass 
by its peak value M M at each field B and the temperature 
by Tm, we observe that all the curves have to merge into 
single one, given by eqs. (5) and (6) thus demonstrating 
a scaling behavior. 

To verify cq. (4), we use measurements of xac(T) in 
CeRu2Si2 at magnetic field B = 0.02 mT at which this 
HF metal demonstrates the NFL behavior [22]. It is seen 
from fig. 2 that eq. (4) gives good description of the facts 
in the extremely wide range of temperatures. The inset 
of fig. 2 exhibits a fit for M^(y) extracted from mea- 
surements of xac(T) at different magnetic fields, clearly 
indicating that the function given by eq. (5) represents a 
good approximation for M^(y). 
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Fig. 2: Temperature dependence of the ac susceptibility xac 
for CeRu2Si2- The solid curve is a fit for the data shown by 
the triangles at B = 0.02 mT and represented by the func- 
tion x{T) = a/VT given by eq. (4) with a being a fitting 
parameter. Inset shows the normalized effective mass versus 
normalized temperature y extracted from xac measured at 
different fields as indicated in the inset [22]. The solid curve 
traces the universal behavior of M^r(y) determined by eq. (5). 
Parameters ci and C2 are adjusted to fit the average behavior 
of the normalized effective mass M^(y). 

Mft{y) extracted from the entropy S(T)/T and magne- 
tization M measurements on the 3 He film [23] at different 
densities x < xfc is reported in the left panel of fig. 3. 
In the same panel, the data extracted from the heat ca- 
pacity of the ferromagnet CePdo.2Rh .s [24] and the AC 
magnetic susceptibility of the paramagnet CcRu2Si2 [22] 
are plotted for different magnetic fields. It is seen that the 
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universal behavior of the normalized effective mass given 
by eq. (5) and shown by the solid curve is in accord with 
the experimental facts. All 2D 3 He substances are located 
at x < xpc, where the system progressively disrupts its 
LFL behavior at elevated temperatures. In that case the 
control parameter, driving the system towards its criti- 
cal point Xpc is merely the number density x. It is seen 
that the behavior of M^(y), extracted from S{T)/T and 
magnetization M of 2D 3 He looks very much like that of 
3D HF compounds. In the right panel of fig. 3, the nor- 
malized data on C(y), S(y), yx{y) and M(y) + yx(y) ex- 
tracted from data collected on CePdi_ x Rh x [24] , 3 He [23], 
CeRu2Si2 [22], and YbRu2Si2 [11] respectively are pre- 
sented. Note that in the case of YbRu2Si2, the variable 
y = B^b/Tm. As seen from eq. (5), this representation 
of the variable y is correct, and Bps makes sense of the 
variable, while the temperature is a fixed parameter. All 
the data show a kink at y > 1 taking place as soon as the 
system enters the transition region from the LFL state to 
the NFL one. Again, we conclude that the presence of the 
kink is mainly determined by the behavior of the effective 
mass at the transition region rather then by the critical 
fluctuations or Kondo scales [11]. 
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Fig. 3: The left panel. The normalized effective mass ver- 
sus the normalized temperature y = T/Tm- The behavior of 
M"^ is extracted from measurements of S(T) /T and magnetiza- 
tion M on 2D 3 He [23]), from ac susceptibility xac(T) collected 
on CeRu 2 Si 2 [22] and from C{T)/T collected on CePdi_ x Rh x 
[24]. The data are collected at different densities and mag- 
netic fields as shown in the left bottom corner. The solid curve 
traces the universal behavior of the normalized effective mass 
determined by eq. (5). Parameters Ci and c 2 are adjusted for 
\n{T n ,B) at B = 0.94 mT. The right panel. The normal- 
ized specific heat C{y) of CePdi_ x Rh x at different magnetic 
fields B, normalized entropy S(y) of 3 He at different number 
densities x, and the normalized yx(y) at B = 0.94 mT versus 
normalized temperature y are shown. The upright triangles de- 
pict the normalized 'average' magnetization M + B\ collected 
on YbRu 2 Si 2 [11]. The kink in all the data is clearly seen in 
the transition region y > 1. The solid curve represents yM^(y) 
with parameters ci and c 2 adjusted for magnetic susceptibility 
of CeRu 2 Si 2 at B = 0.94 mT. 



By definition, MR is given by 
P(B,T) 



p mr (B,T) 



p{0,T) 



(9) 



We apply eq. (9) to study MR of strongly correlated elec- 
tron liquid versus temperature T as a function of magnetic 
field B, The resistivity p(B, T) is 

p(B,T) = p + Ap(B,T) + Ap L (B,T), (10) 

where Ap(B,T) = c(M*(B,T)) 2 T 2 , and c is a constant, 
and the classical contribution A/9^(T, B) to MR due to 
orbital motion of carriers induced by the Lorentz force 
obeys the Kohler's rule [3]: 



Ap L {B,T) / fx B B 
p(0,T) \ P (0,T) 



(11) 



Function F is determined by the details of metal. We 
note that Ap^B) <C p(0, T) as it is assumed in the weak- 
field approximation. Suppose that the temperature is not 
very low, so that p < Ap(B = 0, T), and B > B c0 . 
Substituting (10) in (9), we find that 
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Consider the qualitative behavior of MR described by 
eq. (12) as a function of B at a certain temperature 
T = Tq. In weak magnetic fields, when To > Ti/ 2 and the 
system exhibits NFL regime (see fig. 1), the main con- 
tribution to MR is made by the term Ap L (B), because 
the effective mass is independent of the applied magnetic 
field. Hence, \M*(B,T) - M*(0,T)\/M*(0,T) < 1 and 
the leading contribution is made by Apl(B). As a re- 
sult, MR is an increasing function of B. When B be- 
comes so high that Tm{B) ~ psB ~ To, the difference 
(M*(B,T) - M*(0,T)) becomes negative and MR as a 
function of B reaches its maximum value at Tm(B) ~ To- 
As B increases still further, when Tm{B) > Tq, the effec- 
tive mass M*(B,T) becomes a decreasing function of the 
magnetic field, as follows from eq. (7). As B increases, 
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and the magnetoresistance, being a decreasing function of 
B, is negative. 

Now study the behavior of MR as a function of T at a 
certain value Bq of magnetic field. At low temperatures 
T < Tm(B ), it follows from eqs. (5) and (7) that 



M*(B ) 
M*(T) 



« 1, 



(14) 



and it follows from eq. (13) that p mr (J5o,T) ~ —1, be- 
cause Ap L (B)/p(0,T) <C 1. We note that B must be rel- 
atively high to guarantee that M*(B )/M*(T) < 1. As 



p-4 



Magnetoresistance of the heavy-fermion metal CeCoIn 5 



the temperature increases, MR increases, remaining nega- 
tive. At T ~ Tjvf(-Bo), MR is approximately zero, because 
M*{B Q ) ~ M*(T) and p(B ,T) ~ p(0,T) at this point. 
This allows us to conclude that the change of the temper- 
ature dependence of resistivity p(Bo, T) from quadratic to 
linear manifests itself in the transition from negative MR 
to positive. One can also say that the transition takes 
place when the system transits from the LFL behavior to 
the NFL one. At T > Tm(Bq), the leading contribution 
to MR is made by Apl(Bq) and MR reaches its maxi- 
mum. At Tm(Bq) -C T, MR is a decreasing function of 
the temperature, because 



\M*(B,T) - M*(0,T)\ 
M*(0,T) 



«1, 



(15) 



and p mr (B ,T) < 1. 

The both transitions [25] (from positive MR to negative 
MR with increasing B at a fixed temperature T and from 
negative MR to positive MR with increasing T at a fixed 
value of B) have been detected in measurements of the 
resistivity of CeCoIn 5 in a magnetic field [1]. 

Let us turn to quantitative analysis of MR. As it was 
mentioned above, we can safely assume that the classi- 
cal contribution Api(B,T) to MR is small as compared 
with Ap(B,T). This fact allows us to make our analysis 
and results transparent and simple since the behavior of 
A/9l(_B ) is not known in the case of HF metals. Consider 
the ratio R p = p(B,T)/ p(0,T) and assume for a while 
that the residual resistance po is small in comparison with 
the temperature dependent terms. Taking into account 
cq. (10) and that p(0,T) cx T, we obtain from eq. (12) 
that 



R P = Pr, 



1 = 



P(B,T) 



oc T{M*(B,T)f 



(16) 



It follows from eqs. (5) and (16) that the ratio R p reaches 
its maximum value R P M at some temperature Ta m ~ Tm- 
If the ratio is measured in terms of its maximum value 
R P M and T is measured in terms of Tft m ~ Tm then it is 
seen from eqs. (5), (6) and (16) that the normalized ratio 
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becomes a universal function of the only variable y = 
T/Tjtm. To verify eq. (17), we use MR obtained in mea- 
surements on CcCoIn 5 , see fig. 1(b) of Ref. [1]. The results 
of the normalization procedure of MR are reported in fig. 
4. It is clearly seen that the data collapse into the same 
curve, indicating that MR well obeys the scaling behavior 
given by cq. (17). This scaling behavior obtained directly 
from the experimental facts is a vivid proof that MR is 
predominantly governed by the effective mass M*(B,T). 

Now we are in position to calculate R p v (y) given by eq. 
(17). Using eq. (5) to parametrize M^(y), we extract pa- 
rameters Ci and C2 from measurements of the magnetic ac 
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Fig. 4: The normalized ratio R p N {y) given by eq. (17) ver- 
sus normalized temperature y — T/TR m . The normalized ra- 
tio were extracted from MR shown in fig. 6 and collected on 
CeCoIns at fixed magnetic fields B [1] listed in the right upper 
corner. The starred line represents our calculations based on 
eqs. (5) and (17) with the parameters extracted from ac sus- 
ceptibility of CeRu2Si2 (see the caption to fig. 2). The solid 
line displays our calculations based on eqs. (6) and (17); only 
one parameter was used to fit the data, while the other were 
extracted from the ac susceptibility measured on CeRu2Si2. 



susceptibility \ on CeRu2Si2 [22] and apply eq. (17) to cal- 
culate the normalized ratio. It is seen that the calculations 
shown by the starred line in fig. 4 start to deviate from the 
experimental facts at elevated temperatures. To improve 
the description, we employ eq. (6) which describes the 
behavior of the effective mass at elevated temperatures in 
accord with eq. (4) and ensures that at these tempera- 
tures the resistance behaves as p(T) oc T. In fig. 4. the fit 
of R p N {y) by eq. (6) is shown by the solid line. Constant 
C3 is taken as a fitting parameter, while the other were 
extracted from ac susceptibility of CeRu2Si2 as described 
in the caption to fig. 2. 

Before discussing the magnetoresistance p mr (B,T) 
given by eq. (9), we consider the magnetic field depen- 
dencies of both the peak value R m ax(B) of MR and peak 
temperature TR m (J5) at which R max (B) takes place. It 
is possible to use cq. (16) which relates the position and 
value of the peak with the function M*(B, T). To do this, 
we have to take into account the classical contribution 
Apl(B,T) to MR and the residual resistance po which 
prevent Tn m (B) from vanishing and makes R m ax(B) fi- 
nite at B — > B c q. Therefore, MR is a continuous function 
at the quantum critical point B c q in contrast to M*(B, T) 
which peak value diverges and the peak temperature tends 
to zero at the point as it follows from eqs. (7) and (8). 
As a result, we have to substitute B c for B c o and take B c 
as a parameter. Upon modifying cq. (16) by taking into 
account Apl(B,T) and po, we obtain 



T-Rm(B) ~ bi(B — B c ), 



(18) 
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(19) 



b 3 (B - B e )-i + 1 ' 

Here 6 1; 6 2j ^3 and S c are fitting parameters. It is per- 
tinent to note that when deriving eq. (19), eq. (18) was 
employed in substituting (B — B c ) for T. Then, eqs. (18) 
and (19) are not valid at B < B c q when the HF metal ob- 
tains both the antiferromagnetic order and LFL behavior. 
In fig. 5, we show the field dependence of both TR m and 
Umax, extracted from measurements of MR [1]. It is seen 
that both TR, m and i? ma x are well described by eqs. (18) 
and (19) with B c =3.8 T. We note that this value of B c 
is in good agreement with observations obtained from the 
B — T phase diagram of CeCoIn 5 , see fig. 3 of Rcf. [1]. 
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Fig. 5: The peak temperatures TR m (squares) and the peak 
values i? m ax (triangles) versus magnetic field B extracted from 
measurements of MR [1]. The solid lines represent our calcu- 
lations based on eqs. (18) and (19). 

To calculate MR p mr {B,T), we apply eqs. (17) to de- 
scribe the universal behavior of MR, eq. (5) to describe 
the behavior of the effective mass and eqs. (18) and (19) to 
assign the absolute values to MR. Figure 6 shows the cal- 
culated MR versus temperature as a function of magnetic 
field B, together with the facts taken from Rcf. [1]. We 
recall that the contributions coming from Api,(B, T) and 
Po were omitted. As seen from fig. 6, our description of 
the facts is quite good and we conclude that main contri- 
bution to MR comes from the dependence of the effective 
mass on the applied magnetic magnetic field B. 

In summary, we have performed a study of MR of the 
HF metal CeCoIns within the framework of the fermion 
condensation quantum phase transition. Obtained results 
are in good agreement with facts and have allowed us to 
reveal new scaling behavior of MR. 
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